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A theoretical analysis about the influence of the Hertz elastic contact on a three body chain 
collision is presented. In spite of the elastic character of the collision, the final velocity of each 
particles depends on the particular interaction between them. A system involving two elastic spheres 
falling together one in top of the other under the action of gravity, and colliding with an horizontal 
hard wall is studied in detail. The effect of the Hertz contact interaction can be easily put in evidence 
for some particular situations. 



PACS numbers: 05.45.-a, 46.05. +b 



I. INTRODUCTION 



When two balls are dropped together, both vertically 
aligned and hitting the ground at the same time, an in- 
teresting rebound effect is obtained. If the smaller ball 
is above, it may rise surprisingly high into the air. This 
effect may still be enhanced if one drops a chain of balls, 
ordered by their relative weight, the lightest ball at the 
top. 

This system was already discussed in pedagogical arti- 
cles more than 30 years ago, assuming perfect and inde- 
pendent elastic collisions between each pair of balls [1 5]. 
It was considered an attractive and ideal system to study 
elastic collisions. However, if the balls are dropped to- 
gether, all collisions will happen at the same time. Be- 
cause they are not instantaneous, the final velocity dis- 
tribution of the balls will be dependent on the particu- 
lar interactions between them, even in those cases where 
the interactions are elastic, i.e., with no loss of energy. 
In their beautiful article [2], the class of W. G. Harter 
considered at some point non-instantaneous collisions. 
To study them, they introduced a simple quasi-linear 
model, which could be valid if we think not of spheres 
but of uni-dimensional springs. The dynamical differ- 
ential equations were solved using an analog computer 
program. However, they concluded that the quasi-linear 
model did not agree with their experimental results, and 
they proposed a phenomenological non linear model. The 
interaction field was then obtained directly from empir- 
ical data, by dropping a ball onto a painted flat metal 
surface from varying heights. In this way the interac- 
tion could be plotted as a function of the spot size, and 
eventually as a function of the depression of the ball. Al- 
though the strategy to obtain the non linear interaction 
was remarkable, the discussion that followed was not suf- 
ficiently developed. In fact, from the theoretical point of 
view and as a first approximation, the analysis of the 
non-instantaneous multiple collisions effect should take 
into account the interaction between two elastic spheres 
obeying Hooke's law, which is always valid for sufficiently 



small deformations. This problem is known in the liter- 
ature as the Hertz contact [6]. 

In this article, a system involving multiple collisions 
between two different elastic spheres, one on top of the 
other, with the same initial velocity, hitting a planar hard 
horizontal wall is considered. In Section II, the results 
obtained for instantaneous elastic collisions are briefly 
reviewed. Next, the Hertz contact is introduced in the 
mathematical model of the system. At the end, the re- 
sults of this theoretical analysis are presented. The main 
differences between independent and non-instantaneous 
collisions are outlined and discussed. 



II. INDEPENDENT COLLISIONS 

Consider a system with two elastic spheres, one on top 
of the other, falling with the same velocity v, colliding 
with the ground, as represented in Fig. 1. 




FIG. 1. Two different elastic spheres falling with the same 
velocity colliding with the ground. 

Assume a very small gap between the spheres, such 
that all collisions can be treated individually. First, the 
lower sphere hits the rigid wall and returns with the same 
velocity v. Next, the spheres collide elastically. The final 
velocities u\ and ui (of lower and upper sphere respec- 
tively) can be calculated through momentum and energy 
conservation principle, and are independent of the par- 
ticular interaction. If m = mijm,\ is the spheres mass 
ratio, we can write: 
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where the positive sign corresponds to the upward orien- 
tation. The upper sphere may rebound with a velocity 
up to three times greater than its initial one, for a van- 
ishingly small mass ratio. That is, if its mass is much 
smaller than the mass of the lower sphere. 

In fact, the rebound may be enhanced if we add to this 
system extra spheres, all falling with the same initial ve- 
locity v. To study this effect, consider first the elastic 
collision between two spheres with a mass ratio m but 
with different initial velocities. The upper sphere still 
moves with a downward velocity v. The lower sphere has 
now previously collide with another sphere and moves 
upward with a velocity av (a > 1). The final velocities 
are now given by: 
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Naturally, equations (1) and (2) are recovered if a = 1. 
For a system with n spheres, the upper one may achieve 
a final velocity of (2™ — l)v, if all mass ratios are vanish- 
ingly small. 

Another interesting limit to this system, for a certain 
distribution of mass ratios, concerns the case for which 
all lower spheres stop, and only the upper one rises into 
the air. In this situation, all initial energy was transmit- 
ted to only one sphere. In the simple system composed 
of only two spheres, the lower one has null final velocity 
for m = 1/3, whereas the upper sphere rebounds to reach 
twice its initial velocity. 



III. NON INSTANTANEOUS COLLISIONS 

However, if the spheres fall together, the assumption 
of independent collisions may no longer be accurate. For 
the system represented in Fig. 1, both collisions - of the 
lower sphere with the ground and of the two spheres - will 
happen at the same time. To obtain the final velocities, 
after the spheres separation, it is important to know the 
particular interaction between them. In the following 
subsections, one presents the simplest possible interac- 
tion: it will be assumed that all solids respect Hooke's 
law of elasticity, valid for sufficiently small deformations. 
The dynamical equations for this model will then be es- 
tablished and solved. 



spheres, with radius of curvature respectively of Ri and 
i?2 (see Fig. 2), the elastic potential energy E p depends 
on the combined deformation h as: 
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where the reduced radius r — i?ii? 2 /(i?i + R2) and the 
reduced elastic constant e is dependent on the Young 
moduli Ei, E 2 and on the Poisson coefficients ai, a 2 . It 
is given by: 
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The Young modulus increases with the solid rigid- 
ity, whereas the Poisson coefficient is typically slightly 
smaller than 1/2. 




FIG. 2. Elastic contact between two spheres. 

The elastic energy (5) is valid for a static deformation. 
Nevertheless, it can be considered a good approximation 
for the situation studied in this article if the velocity v 
is much smaller than the sound velocities of the solids 
involved. It is interesting to notice that this approxima- 
tion fits well the experimental potential energy obtained 
by the class of W. G. Harter [2], in the first regime of 
small ball depressions. For large depressions, the energy 
increased with a greater power of h, which is in agree- 
ment with the fact that Hooke's clastic regime only take 
into account the first terms of the elastic energy power 
expansion on the deformation. 

It is instructive to calculate the maximum height of de- 
formation and the contact time of two colliding spheres. 
Suppose v is the relative velocity of the spheres before 
the collision, and /1 = mimil(rri\ + U12) their reduced 
mass. Then, the energy conservation principle reads: 
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The maximum height of deformation occurs when the 
relative velocity is zero and can be written: 



A. The Hertz contact 
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The theory of the elastic contact between solids was 
first studied by H. Hertz, and it can be followed in detail 
in Ref. [6]. For two isotropic and homogeneous elastic 



The collision time corresponds to the time in which the 
deformation goes from to Km and back to again. It 
is given by: 
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These are the relevant length and time for elastic col- 
lision problems, and they will be useful to write down 
dimensionless dynamical equations. 



B. Dynamical equations 

For the system represented in Fig. 1, the total potential 
energy is given by the sum: 

2 1 i 2 1 £ 

Ep = -^emr^h^ + -e 12 rf 2 hf 2 , (10) 

where the first term takes into account the interaction 
between the first sphere and the rigid ground, with in- 
finite Young modulus and radius of curvature, and the 
second term is the interaction between both spheres. The 
deformations depend on the positions X\ and x 2 of the 
centers of the spheres (the ground defining the reference 
position): 

h 01 = H{R 1 - Xl ), (11) 
h 12 =H((R 1 +R 2 ) - {x 2 -xi)). (12) 

The function H{x) = x if x > and H{x) = otherwise. 
The equations of motion m,ij = —dE p /dxi, (i = 1,2), 
will simply be written: 

mixi = eoir^h^ - e 12 r? 2 hl 2 , (13) 
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m 2 x 2 = e 12 rf 2 hl 2 . (14) 
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If one introduces natural units, L — (mit> 2 / e oi?*oi) b and 
T = L/v, which are associated with the maximum height 
of depression and the time of contact between the ground 
and the first sphere, it is possible to define dimensionless 
variables. The equations of motion may be written in a 
simpler form: 

3 3 

xi = - khl 2 , (15) 

3 

mx 2 = khl 2 . (16) 

They only depend on two parameters, 

m = m 2 /m\, (17) 
i i 

k = e 12 rl 2 /e 01 r^. (18) 

Note that k is defined as a ratio of the reduced elastic 
constants and radius. If the infinite rigidity of the pla- 
nar wall is taken into account, we have eoi > e 12 and 
foi > r\ 2 . For this system, one concludes that k < 1. 



C. Numerical resolution 

The set of differential dynamical equations (15) and 
(16) was solved numerically using a simple and stan- 
dard Eulcr method, for different values of the parameters. 
Fig. 3 shows several solutions found for the positions x\ 
and x 2 as functions of time, starting when the two spheres 
initially in contact hit the ground with the same veloc- 
ity v, until they separate with constant final velocities u\ 
and«2. The axis have arbitrary units. The solutions were 
calculated for spheres with the same Young modulus and 
density, but four different mass ratios m = 0.01, 0.5, 1,5 
(wich were obtained choosing appropriate sphere radii). 




m=0.01 t m=0.5 t 




m=l t m=5 t 



FIG. 3. The positions xi and X2 as a function of time, for 
spheres with the same Young modulus and density, and four 
different mass ratios m — 0.01,0.5,1,5. The axis have ar- 
bitrary units. The bold lines represent the ground and the 
center of the spheres. In thinner lines, x\ ± R\ and X2 ± R2 
are also represented, to indicate in which periods of time the 
interaction exists. 

Surprisingly, for a very small mass ratio, the spheres 
have almost the same velocity u\ ~ u 2 w v after the 
rebound. This situation contrasts with the results ob- 
tained considering independent collisions, for which the 
maximum velocity gain would be expected. In fact, if the 
mass ratio m is negligible, it is possible to see from Eq. 16 
that the deformation h\ 2 will also be very small, exhert- 
ing almost no influence in the motion of the lower sphere 
(see Eq. 15). The latter sphere rebounds with the wall 
and returns with velocity u\ w v. Since h\ 2 w R 2 w 0, 
it may be concluded that R\ ~ x 2 — x\. Then x,\ ~ x 2 , 
which means that both spheres stick together during the 
collision as if they were one whole body [7] . 

It is also interesting to notice that for mass ratios 
greater than one (in the case of the figure with m = 5) 
the three body collision gets more complex, with the bot- 
tom sphere making several rebounds between the ground 
and the other sphere. 

It is important to mention that the final velocities do 
not depend on the amplitude of the spheres' deformations 
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during the collisions. If both spheres were more rigid, but 
keeping the same ratio k, the deformations could be very 
small, and still the final velocities would be equal to the 
ones represented in Fig. 3. 

In Fig. 4, the bold lines represent the final velocity dis- 
tributions u\ and u 2 as functions of to < 1, for different 
values of E 2 jE\ = 0.1, 0.5, 1, 10. The thin lines show the 
results of equations (1) and (2), considering independent 
collisions. 




FIG. 4. Velocity distributions in and U2 as functions of 
m < 1. The thin lines correspond to the results obtained for 
independent collisions and they are independent of the ratio 
between the Young moduli. In bold, the velocities obtained 
considering the Hertz contact for E 2 /Ex = 0.1, 0.5, 1, 10. 

The independent collisions results are recovered for 
small values of the ratio E 2 /E\, and consequently for 
small k. This is particularly true for the final velocity 
distribution of the lower sphere u\. In fact, if k — > 0, 
it is possible to see from Eq. 15 that in a first moment, 
the motion of the lower sphere will depend mainly on 
its interaction with the wall - the deformation h\ 2 be- 
tween both spheres will have a negligible contribution 
when it is compared with the deformation /ioi between 
the lower sphere and the wall. However, if the collision 
time between both spheres is large enough (as it is shown 
in Eq. 9, the collision time increases for small k), after 
the lower sphere-wall interaction, the terms containing 
the deformation h\ 2 will be the only ones present in the 
equations of motion. Therefore, the collisions can be con- 
sidered as almost independent and subsequential. Never- 
theless, as the final velocity u 2 is concerned, the conver- 
gence does not occur for all values of to: it is possible to 
see in Fig. 4 (E 2 jE\ = 0.1) that the final speed u 2 — > v 
as the mass ratio to — > 0. 

As E 2 /Ei (or k) increases, small changes can be seen 
in the final velocity distributions, and the independent 
collisions approximation is no longer appliable. These 
changes arc limited because the infinite rigidity assump- 
tion requires k < 1. In fact, the velocity distributions for 
E 2 jE\ >• 1 is not very different from the one represented 



for E 2 /E 1 = 10 

The left part of Fig. 5 represent the rebound maximum 
speed u 2m as a function of k < 1. On the right, the mass 
ratio rriM corresponding to the maximum rebound speed 
is also plotted, tom decreases rapidly to zero as k — > 0. 
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FIG. 5. Maximum rebound speed U2 M and corresponding 
mass ratio ttim as a function of k < 1 (according to the infinite 
rigidity of the wall assumption). 

The discontinuity on the slope observed for the final ve- 
locity Mi in Fig. 4 (E 2 jE\ = 0.1) corresponds to the criti- 
cal value m c for which the lower sphere stops completely. 
If to > to c , this sphere rebounds twice off the ground, be- 
fore it reaches the final velocity represented in the figure. 
If the ratio between the Young moduli increases, the dis- 
continuity changes first its position {E 2 /E\ — 0.5), and 
eventually it disappears [E 2 jE\ = 1 and larger values). 

However, these discontinuities will happen again for 
larger values of to, as the lower sphere rebounds more 
and more in between the ground and the upper sphere. 
This curious feature can be seen in Fig. 6, which repre- 
sents the velocity distributions u\ and u 2 as functions of 
1/to < 1, for E 2 /E\ = 1. As 1/to approaches 0, i.e., 
when the lower sphere becomes much lighter than the 
top one, the number of collisions increases to infinity. 




FIG. 6. Velocity distributions ui and 112 as functions of 
1/m < 1, for E2/E1 — 1. These results correspond to the re- 
sults obtained only considering non-instantaneous collisions. 

IV. DISCUSSION AND CONCLUSIONS 

In this article, the influence of the Hertz contact on 
multiple chain collisions was studied in detail. It was 
shown that the independent collisions approximation is 
no longer valid for small values of the mass ratio to and 
large values of the rigidity ratio k. 

To understand the results obtained here, it is impor- 
tant to consider the typical two-body times of collision 
(see Eq. 9). First, suppose to is of the order of unity. If 
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k is small, the two-body collision time between the lower 
sphere and the wall, tqi, is smaller than the two-body 
collision time between both spheres, t\i- As k — ► 0, the 
first collision occurs almost instantaneously, and subse- 
quently, the two spheres collide. Thus, the results ob- 
tained considering this limit approach the independent 
collisions results. 

However, as k increases, tq± also increases. The col- 
lisions become more evolved, and it is difficult to treat 
them separately. 

Despite its complexity, the final velocity distributions 
Mi and U2 approach each other for large k and small to. 
In this limit (m — > 0), as it was prooved through the 
analysis of the equations of motion, both spheres behave 
as if they belonged to the same body, changing their ve- 
locities together as they interact with the ground. At 
the end, both spheres rebound with equal final velocities 
Ul « U2 ~ v. 

The results obtained with this model were summarized 
in Fig. 5, which clearly indicates the best parameters m 
and k to experimentally put in evidence the influence of 
the Hertz contact on chain elastic collisions. 
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